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ABSTRACT 

We measure the matter probability distribution function (PDF) via counts in cells in a volume 
limited subsample of the Sloan Digital Sky Survey Luminous Red Galaxy Catalog on scales from 
30 /i -1 Mpc to 150 h~ 1 Mpc and estimate the linear Integrated Sachs- Wolfe effect produced by super- 
voids and superclusters in the tail of the PDF. 

We characterize the PDF by the variance, S3, and S4, and study in simulations the systematic 
effects due to finite volume, survey shape and redshift distortion. We compare our measurement to 
the prediction of ACDM with linear bias and find a good agreement . 

We use the moments to approximate the tail of the PDF with analytic fu nctions. A sim ple Gaussian 
model for the superstructures appears to be consistent with the claim by IGranett et al\ that density 
fluctuations on 100 /i _1 Mpc scales produce hot and cold spots with AT s» lOfiK on the cosmic 
microwave background. 

Subject headings: cosmic microwave background — large-scale structure of universe — methods: sta- 
tistical 



1. INTRODUCTION 

After last scattering photons traveled through mostly 
neutral media. Although radiation and matter are not 
strongly coupled, there is still a secondary signal due to 
large scale structure on top of the primary fluctuations of 
the cosmic microwave background (CMB ) radiation. The 
Integ rated Sachs- Wolfe (ISW) effect (jSachs and Wolfd 
Il967t ) accounts for most of the secondary anisotropies 
for low multipoles (|Hu and Dodelsor] 120021 ) . As the ex- 
pansion of the universe accelerates, gravitational poten- 
tial wells and hills decay. Photons traversing these get 
blueshifted or redshifted. 

Due to its weak signal, ISW detection is very challeng- 
ing. Cross-correlating galaxy surveys w ith CMB maps 
yield r esults from ma r ginall y significant dScranton et all 
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Raccanelli et al. 2008 
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Sawangwit et al. 2010) to 4.5a de 



tections ( Giannanto nio et 



2008). The 



higher significance was achieved by a joint analysis of 
surveys. Other techniques focusing on the signal from 
discr ete objects can reach up to 4.5g from a single sur- 
vey (jMcEwen et al\ 120081: IGranett et a?.ll2008D . 

The ISW effect can be a unique probe of dark en- 
ergy if well-measured. From cross-correlation measure- 
ments and the Wilkinson Microwave Anisotro py Probe 
( WMAP) power spectrum ([Bennett et a/Jl2003h it has al- 
ready been shown t hat it is possible to constrain cosmo- 
logica l parameters (|Giannantonio et all l2008t IHo et al] 
2008). Despite the fact that the detection of the signal 
from discrete objects has higher significance, they can- 
not be used for parameter estimation due to the lack of 
simple quantitative models. 

Further motivation for studying super structures stems 
from anomalies in the low 1 mode s of the CMB 
(jTegmark et q/J|200l ICopi et a/.ll200l . Ilnoue and Silkl 
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(|2007|) calculate the effect of large, dust filled, compen- 
sated voids in the local universe. They were successful in 
explaining the observed CMB anomalies but these voids, 
due to their size and depth, do not fit into the widely 
accepted picture of clustering. They assume extra power 
on large scales in the matter power spectrum. When 
subtr acting the estimated l ocal I SW signal from CMB 
maps, iFrancis and P eacock (2010) found that the signif- 
icance of the anomalies decreased. 

In this paper, our principal goal is to give an estimate 
of the ISW signal coming from large overdense or under- 
dense regions (superclusters or supervoids). To achieve 
this we have to deal with two separate problems: the 
density of the extreme fluctuations and the ISW effect 
associated with them. In Section^ we present our mea- 
surement of counts in cells (CIC) in the Sloan Digital 
Sky Server (SDSS) Luminous Red Galaxy (LRG) Cata- 
log. From CIC we derive the first few moments (the vari- 
ance, 5*3, and S4) of the matter probability distribution 
function (PDF) and compare them to their theoretical 
values. In Section [3j we use these to estimate the PDF 
focusing on its tail, since an enhanced tail could explain a 
strong ISW signal from super structures. We use a simple 
Gaussian model to derive an expression for their profile 
and the potential. We com pare their estimated I SW sig- 
nal to the elusive results of IGranett et al\ (|2008f ). More 
about their ISW measurement can be found in Section 
[3] of this paper. In Section |4j we summarize and discuss 
our results. 

2. MEASUREMENTS 

The hydrodynamical model of the universe is based 
on the assumption that the observed galaxy distribution 
is a Poisson-sampled version of a continuous field. Fur- 
thermore, this conti nuous field is a realization of a ran- 
dom field (see, e.g.. lPeebleslll980l ). Measuring CIC is a 
well-established method to estimate its P DF (see, e.g., 
iColombi et ~al\ l2000t Isiapudi et al\ [2000). In this Sec- 
tion first we give a brief summary of this method and 
description of the data and the algorithm we use in our 
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analysis. Then we compare our findings to the predic- 
tions of ACDM. 

2.1. CIC and the PDF 

In the following we deal with two random fields, one 
corresponding to dark matter, and the other its biased 
version, the galaxy field. In our notation the terms that 
refer to these are " matter PDF" , " matter field" , " galaxy 
PDF" and "galaxy field". 

The matter PDF for cel ls is fully give n by its cumu- 
lants, possibly normalized (|Peebleslll980h : 

?n = V^Iv-Sv (S(xi)...S(x n )) c dx 1 ...dx n , (1) 

S n = jkr, (2) 

42 

where V is the volume of a cell and the subscript c refers 
to connected moments. When £ 2 << 1, S n 's are around 
unity according to perturbation theory. Thus, S3, the 
skewness, and S4, the kurtosis represent the lowest or- 
der correction to a Gaussian distribution. Deriving these 
quantities from galaxy counts can be done in two steps: 
first, by using factorial moments to get the moments of 
the u nderlying galaxy PDF (see e.g. iSzapudi and Szalavl 
|1993]), second, by using some galaxy-dark matter bias- 
ing scheme to transform the galaxy PDF into the matter 
PDF. 

Estimators for the variance (£2)1 ^3 an d Si are: 

- _ (n*)-(n) 

{N y - 1, (3) 

g 3 = (F 3 -3F 2 F 1 +2F 1 3 )/F 1 3 ,^ 

_ (F i -4F 3 F 1 + 12F 2 Fj- 6F t 4 - 3F 2 2 ) / F* ,^ 

where 

F n = (N(N-l)...(N-n + l)). (6) 

Ergodicity ensures that these ensemble averages can be 
calculated from a single, ideally large, volume limited 
survey. The caveats, arising when we depart from the 
ideal case, are discussed later in this Section. 

In this paper we use cosmol ogical parameters ta ken at 
their best-fit WMAP values (jSpergel et aLll2007f ). For 
the bias we fit the simplest, deterministic, local, linear 
model: 

5 g = bS. (7) 

This is generally a good approximation on quasi-linear 
scales. Its validity is tested in Section l2~5l 

2.2. The Data 

Among the spectroscopic galaxy surveys available to- 
day the Seventh Data Release of the Sloan Digital 
Sky Survey (SPSS D R7) covers the largest volume 
(jAbazajian et all 120091 ) . The LRG sample is generally 
regarded as a good cosmol ogical probe. P r operti es of 
the LRGs can be found in lEisenstein et al\ (|2001l ). In 
practice, flags in the Sloan database identify these galax- 
ies. For CIC one needs a volume limited sample which 
can be obtained by magnitude and redshift cuts. We re- 
stricted our LRG sample to redshifts between 0.24 and 



0.31 and fc-corrected absolute magnitudes between —22.3 
and —24.3 in the r band. We used values from the Photoz 
and SpecObjAll tables on the SkyServer Web site. By ex- 
cluding the three stripes in the Southern Galactic Cap we 
were left with 21613 galaxies. After converting redshift 
into comoving radial distance by using the best WMAP 
cosmological parameters, the angle-averaged density ap- 
pears to be uniform with fluctuations consistent with 
Poisson noise. The selection function of a similar data 
set is plotted in Figure 12 in lEisenstein et a!\ (|2001D . 

2.3. The Algorithm 

From R.A., decl., and z coordinates, we calculated co- 
moving Cartesian coordinates. Then we placed a rect- 
angular grid over the sample. Since in this arrangement 
cubical cells are readily accessible, we chose to measure 
CIC in cubes. The survey mask, however, has a com- 
plex shape, usually given by spherical polygons. Holes 
and the irregular boundary cause unwanted edge effects 
which could bias the results in a complicated way. 

To tackle this problem, first we took a cube-shaped 
region encompassing the survey. Then we created two 
negatives by filling the parts in the mask and outside the 
survey area with dummy galaxies from a Poisson point 
process; one with the average density of the survey and 
one with hundred times that density. We added the first 
negative to the survey to fill the holes. We measured 
CIC in this and in the second negative in parallel. Since 
its density is large, the counts from the second negative 
provide a good measure of the overlap of the cells with 
the survey geometry. We ignored any cell that had more 
than 10% of its volume outside, which corresponded to 
having a galaxy count in the second negative larger than 
100 X average density X volume of cell X 0.1 . In this 
work we used MANGLE (|Swanson et a/.ll2008D to check 
whether an object was inside the mask. 

2.4. The Systematics 

The shape of the survey and diluting the data with 
a Poisson point process introduce systematic bias into 
our measurements of CIC. In order to assess its level 
we studied simulations. We created mock catalogs with 
a second-order Lagrangian (2LPT) code (jCrocce et al\ 
I2006D . We created 100 mock catalogs in 2500 /i _1 Mpc 
cubes, then we used these to create another set of mocks 
by applying the mask of the spectroscopic survey. The 
galaxies were downsampled in every case to match the 
average density with that of the data. With these and 
the two negatives described in Section [2~31 we were able 
to measure CIC in three different arrangements. 

(i) Ideal, large, cubic-shaped simulations. 

(ii) Variant of (0). We only took into account a cell 
when at least 90% of its volume lay inside the survey 
area. 

(iii) Dropping galaxies outside the survey area and fill- 
ing it with the dummy galaxies of the negative to 
preserve the average density. (The 90% rule still 
applies) 

First we estimated the cosmic bias due to the survey 
shape and volume. This qu estion has been studied exten - 
sively in the past (see, e.g.. ISzapudi and Colombil ll996'). 
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We measured the variance, S3 and S4 in the first two 
arrangements (0 and (H3). In Figure [TJ the ratio, 6A/A, 
is plotted for each of these quantities, where A is the 
quantity measured in arrangement (0) . Error bars were 
estimated from the scatter around the average. The error 
of the average is plotted, so the error of a single measure- 
ment is ten times larger. In the case of the variance this 
is the well-understood integral constraint problem and 
the ratio does not exceed a couple of percents even at 
the largest scale. For S3 the ratio is consistent with 1 
but for S4 the difference from 1 is not negligible even at 
relatively small scales. However, as we show later, this 
bias is still small compared to the cosmic error. 

The bias, caused by the data having been diluted with 
a random sample, can be understood in case that evenly 
distributed holes comprise the mask. The resulting cata- 
log can be considered the linear combination of the Pois- 
son sampled galaxy and a constant field. This constant 
field is the random sample that fills the mask: 

P = Pgalaxy T" Prandom- (8) 

Subsequently the density contrast can be written as 

<5 = ySgalaxy + (1 - y)S rand om (9) 

with 



y 



(10) 



Since S ranc iom is zero the moments of S are proportional 
to the moments of S^ata '■ 

I = V 2 laa. 



galaxy ) 

q „.— n+2 qgalaxy 

On — y O n 



(ii) 

(12) 

We used arrangement (JTTJ) and ([ml) to express y ac- 
cording to Equation (fTTI) an d Equation ([T2]) (Figure HJ. 
The measured values are consistent with our assumption, 
Equation (fTU)) . The robustness of this simple model is 
due to the fact that the variance and S n are insensitive 
to small changes in the cell shape (e.g., in the case of 
a p ower-law corre lat ion function the S„s are constants, 
see IPeeblesI 119801 or IBoschan et all 119941: ISzapudil 119981 
for a study of on this). This bias can be corrected for by 
measuring y directly. 

As the next step we added redshift distortions to the 
mock catalogs. The effect on the variance is expected 
to be similar to the effect on the monopole of the two- 
point fu nction. S3 a nd Sa are affected less according 
to, e.g., [Hivon et al\ (11995 | ). While the Kaiser formula 
(|Kaiserl ll987t IHamiltonl 119921 ) gives a good description 
of this in the linear regime. The lowest order of the 
th ree-point function in Fouri er space has been worked out 
by iScoccimarro et al\ (|1999t ) but Fourier transforming it 
back to redshift space is infeasible in general. Higher 
moments are gradually harder to compute. For these 
reasons we follow a phenomenological approach. In Fig- 
ure [31 the ratios of redshift and real space values of the 
variance, S3 and S4 are plotted. The thick line is the 
predicted amplification from the Kaiser formula for the 
variance: 



with 



^ s = (l + 2//3 + / 2 /5)£, 

f = n°- 6 /b, 



(13) 
(14) 
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Fig. 1. — Bias due to the finite volume of the survey as a function 
of the cell size. 
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Fig. 2. — y as a function of cell size from Equation JTTJ and H12I I. 



where the subscript RS stands for redshift space. In the 
simulations b is 1. In this paper we assume that the effect 
of redshift distortion is small compared to the cosmic 
errors in the case of S3 and S4, and that the variance is 
amplified according to Equation (JT3J) . 

In practice if the systematic bias is small compared to 
the cosmic error then it is negligible. In FigureUJthe total 
systematic bias after corrections according to Equation 
(fTTj) and ([T2"j) and the cosmic error are plotted. It can be 
concluded that the proposed corrections are sufficient to 
measure the variance, S3, and S4 with an error that is not 
significantly different from the cosmic error. This plot 
also tells us that the signal-to-noise ratio drops below 1 
around 100 /i -1 Mpc for S3 and around 50 /i _1 Mpc for 
S4, so they cannot be measured reliably beyond these 
scales. 

2.5. The Variance, S3 and S4 
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Fig. 3. — Effect of rcdshift distortion on the variance (£), S3, and 
S4 as a function of the cell size. 



1.0 
0.5 




40 60 80 100 120 140 




L(Mpc/h) 

Fig. 4. — Total expected syste matic bias is plotted after correc- 
tions as discussed in Section 12.41 The shaded region represents the 
cosmic error. 

We measured CIC in the SDSS DR7 spectroscopic 
LRG sample and compare the prediction of ACDM to 
our results. 

For measuring the variance, S3 and S4 we followed the 
procedure outlined in Section 12.31 and we corrected for 
the systematic bias as given by Equation (fTTj) and (TT2"j) . 
We determined the y parameter from the simulations. 

Additionally, in the case of the real data one has to as- 
sume a galaxy-dark-matter biasing scheme. We used the 
simplest linear local model as given by Equation (JTJ) . We 
computed the bias parameter from fitting the variance 
and testing its consistency on S3 and S4. We defined the 
following chi-square: 

X 2 (b) = (h-U^C-H^-Uib)). (15) 

Here £ d and £ th stand for the measured and the theo- 
retical variance. For the theory we used the real space 
linear model and we assumed that it transforms to red- 



FlG. 5. — Measured variance, £3, a nd S 4 from the SDSS spec- 
troscopic LRG subsample (see Section 12.21 for a description of the 
data) with the theoretical predictions (solid line) vs. the cell size. 
L is the size of a cubic cell. 

shift space as the monopole of the two-point function (see 
Equation (fTU)) ). The measured variance was rescaled as 
in Equation ([lip and extrapolated to present day ( z = 0) 
using the growth function (see, e.g., lDodel son 2003). The 
covariance matrix was calculated from mock catalogs de- 
scribed in Section l2~4l 

Ck = lvb£n(£ -ZKg-Ti), (16) 
^=££X (17) 

where the superscript n refers to the nth simulation and 
N is the total number of simulations, in this case 100. 
These simulations are in redshift space and with bias 
equal to 1. After finding the minimum of the chi-square, 
the covariance matrix was rescaled according to Equation 

For our fit we used the range from 30 h 1 Mpc to 
150 /i _1 Mpc where the linear theory is generally as- 
sumed to be valid. The result is b = 2.14+q'jJ with 
lcr uncertainty , whic h is consistent with findings of 
lOkumura et al\ (|2008D . who used a very similar data set. 
The best fitting variance is plotted with the data on the 
upper panel of Figure [5] When we changed the bound- 
aries of the range to 50 h Mpc and 130 h~ 1 Mpc, we 
found that the change in the best b was consistent, only 
+0.03. 

We also measured S3 and S4 and applied Equa- 
tion (jT2"j) . In Figure [5] these are plotted along with 
the prediction of linear ACDM (|Juszkiewicz et a/]|1993l ; 
lBernardeau1ll994| ): 

bS 3 = f+Ji, (18) 
b 2 S± = ^ + f 7i + hi ~ §72 (19) 

with 

7« = ZT5p- (20) 

The average of the correlation function, £, was calcu- 
lated via Monte Carlo simulations for cells of 30 /i _1 Mpc 



5 





s 




^3 


A Si 




AS 4 


30 


0.980 


0.046 


1.951 


0.238 


7.176 


2.225 


40 


0.574 


0.034 


1.916 


0.273 


6.701 


2.234 


50 


0.370 


0.027 


1.804 


0.343 


5.821 


2.770 


60 


0.256 


0.022 


1.613 


0.452 


4.191 


3.581 


70 


0.187 


0.019 


1.404 


0.586 






80 


0.144 


0.017 


1.179 


0.755 






90 


0.115 


0.015 


0.942 


0.951 






100 


0.094 


0.013 










110 


0.078 


0.011 










120 


0.066 


0.010 










130 


0.056 


0.009 










140 


0.048 


0.008 










150 


0.040 


0.008 











TABLE 1 

The Numerical Values of £, S3, and 54 with la 
Uncertainty. We leave the fields blank when the 
signal/noise ratio is less than one. 
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A 


P 




0.0006 


0.05 




0.13 


0.80 




1. 


0.83 


Joint 


1. 


0.85 



TABLE 2 

Testing the Lin ear Model with Shrinka ge Techique as in 

IPOPE AND SZAPUDJ (120081) . 

+A,40 /i^Mpc +A,...,150 ft^Mpc +A, where A is 
or ±3. The 7 values were estimated using discrete deriva- 
tives of £. The data values from the plot are collected in 
Table Q] 

We tested the goodness of the theory by calculat- 
ing a covariance matrix with the shrinkage technique 
(|Pope and Szapudl 120081 ) with a diagonal of the empiri- 
cal covariance matrix (C) as the target covariance matrix 
(T): 

C i:j = XTij + (1 - X)Cij (21) 

A rec ipe to calculate A is given in iPope and SzapudU 
( 2008). This method ensures that we get a well behaving 
covariance matrix. We estimated the significance for the 
variance, S3 and S4 separately and jointly. The results 
are in Table [2j showing a good overall agreement with 
our ACDM model with linear bias. Visually, the S n 's 
appear to be slightly larger than expected on smaller 
scales, however, this is not statistically significant at all. 
On large scales, there appears to be a slight excess power 
on 130 — 150 /i _1 Mpc scales in the variance (not appar- 
ent in Figure [5]) . This is not large enough to influence 
the ISW effect, and its significance is only 2a according 
to Table [2] In summary, ACDM is a good fit for all the 
moments we measured. Note that the shrinkage estima- 
tor (A = 0.0006) gave identical results for the covariance 
matrix of the variance to that given by Equation (fTTj) . 

3. THE ISW EFFECT OF EXTREME FLUCTUATIONS 

Here, we explore the possibility that the largest fluc- 
tuations in the linear matter density cause detectable 
a nomalies on th e CMB through the linear ISW effect. 

iGranett et all (|200 8) identified large underdense and 
overdense regions (supervoids and superclusters) in the 



SDSS DR4 photometric LRG sample. They stacked im- 
ages, cut out from the CMB, centered on the directions of 
50 supervoids and 50 superclusters found with the high- 
est significance. 

The signal is consistently present in every frequency 
band, so it is likely to have a cosmological origin. In 
this Section we put under scrutiny the tail of the den- 
sity distribution and the density profile of supervoids or 
superclusters. The question we ask is: what is the ex- 
pected ISW signal produced by the 50 objects with the 
highes t and lowest densitie s in a survey similar to the 
one in IGranett et oli (pOOl )? 

In order to answer this we have to have an estimate for 
tail of the matter PDF. We show that by using simple 
analytic functions to approximate the galaxy PDF this 
can be done robustly for our purposes. We also revise the 
way in which the ISW signal from these density extrema 
is estimated. 

3.1. The Matter PDF 

If the galaxy PDF is known, a simple convolution with 
a Poisson-distribution gives the galaxy counts: 

P(N) = J <"> ^ )N e-( N )^P(6 g )d5 g . (22) 

On large scales one can approximate P(S g ) as Gaus- 
sian, lognorm al or second-order Edgeworth expan- 
sion (see, e.g., iKim and Stra uss 1991; iSzapudi and Pan! 
I2004f h The first two depend on the variance only, while 
for the Edgeworth expansion we need S3 and S4 as well: 

P o( S 9) = (23) 
P ^(<y = ^e^+^/^/^/il + 5,) (24) 

M* a )= 7 ^^ /2CT2 (l + £ r^3(W 

+ ^H 4 (S g /a) + ^fH 6 (6 g /a)), (25) 

where a 2 — £, a 2 = y^l + £) and Hi is the ith Hankel- 
function. The above approximations are plotted in Fig- 
ure [5] for galaxy counts in cubes with a linear size of 
100 /i -1 Mpc from the spectroscopic LRG sample. The 
best match around the tails is the Edgeworth expansion, 
while around the mean all of these approximations pro- 
vide qualitatively similar results. The left panel of Figure 
[S] shows the underlying continuous galaxy PDFs, the de- 
convolved CIC distribution, while the right shows them 
convolved with the Poisson-distribution and the galaxy 
counts. 

To estimate the density of the 50 supervoids and 50 
superclusters of a certain size we generated V sur vey /V C eU 
random numbers with the PDFs given above, where 
Vsurvey refers to the volume of the survey and V ce u refers 
to the volume of a superstructure. This approximation 
slightly overestimates the number of independent cells. 
Then we stored the lowest and highest 50 of the num- 
bers. Repeating this several times gives the distribution 
of the extrema. This method, however, is too slow for 
calculating a covariance matrix from hundreds of simu- 
lations. Since we only deal with linear ISW here, as we 
explain it in the next Subsection, we need only the mean 
of the extrema rather than their whole distribution. A 
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Fig. 6. — On the left panel, three analytic functions approximating the galaxy PDF for cells of 1 00 fe~ 1 Mpc are plotted (for details, see 
the text). On the right panel, these are plotted again convolved with a Poisson noise (Equation J22}). The data is the solid line on the 
right panel. 



satisfactory approximation for the mean in the case of 
voids is: 



5 a = P(8 g )8 g dS g 



(26) 



where S max is given by 50V ceU /V survey = J_™" x P(S g )d5 g . 
For clusters the limits of the integration change to J^ . 
In our tests we found the difference to be about a few 
percent. As we use the linear bias model, a further divi- 
sion by the bias yields S. 

3.2. The Profile of Supervoids and Superclusters and 
the ISW Signal 
3.2.1. The Average Profile of Superstructures 

In general the ISW effect is determ ined by an inte- 
gral a long the path of a CMB photon (|Sachs and Wolfe] 
1967): 



AT 2 f , 9$(r(r),r) 

= — tt / dr . 

T c 2 ./ dr 



(27) 



where r denote s the conformal time. 

According to iRudnick et a l. (2007) a simple estimate 
of the linear ISW effect, an undcrdense or an overdense 
spherical region at redshift z causes a 



AT 
T 



M^) 3 (i + 2z)(i + z)- 2 <5 



(28) 



temperature shift at its center on the CMB, where r c is 
the comoving radius of the sphere. In the derivation the 
authors approximated the potential with a top hat which 
implies a compensated void or cluster. Here, instead, we 
propose a profile motivated by Gaussian statistics. For 
a spherically symmetric object the Newtonian potential 



can be calculated easily as 
$(r) 



3n m 

8tt 



Ho 



(29) 



where M(r) = 4ir J Q drf 2 S(f). The density contrast at 
distance r from the center can be obtained with the con- 
dition that the average density inside r c is known: 



P{S in ) 



P(6(r),6. 



• exp 



. exp 



\SC~H 



P(S(r)\8 in ) = P(S(r),S in )/P(S in ) 



(30) 

(31) 
(32) 



where (...) refers to unconditional ensemble averaging 
and Sin is the average density measured inside r c . For 
short we use 8 — (S(r), Si n ) and C = (8®S*) uc . In practice 
(...) can be computed with Monte Carlo simulations 
robustly. From these the expected density is 



(33) 



This profile is not compensated inside a finite radius. 

3.2.2. The Uncertainty of the Profile and the ISW Effect 

It is useful to calculate the uncertainty of the profile in 
order to get an estimate of the uncertainty of the poten- 
tial and the ISW effect. We would like to point out that 
for a correct treatment of the potential one should drop 
the assumption of spherical symmetry. We chose to op- 
timize the accuracy and speed by keeping the spherical 
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Fig. 7. — On the upper panel M(r) = in J Q r drr 2 S(r) is potted 
for a supercluster with Si n = 1 for r c = 62 h~ Mpc. On the lower 
panel we plot the potential (dashed line) and the top hat potential 
(see the text for details). 

approximation. We use the Gaussian model as before: 

P(6(n),6(r 2 ),5 in ) = -y=^ex P ( - ^C^aW) 
P(6(ri),6(r 1 )\S in )=P(5(r 1 ),5(r 1 ),5i n )/P(.Si a ), (35) 

where (5*3 = (S(n), 5(r 2 ), S in ) and C 3 = (63® 63)^. From 
this the covariance between shells at 7*1 and r 2 is: 

Cov(r(l),r(2)) = (S(n)5(r 2 )) - (S( ri ))(S(r 2 )) 



-(5(ri)5 in ) uc (5(r 2 )5 in )^ 



(36) 



On the upper panel of Figure [7] we plot M(r) for r c = 
100/(4 7 r/3) 1/3 h- 1 Mpc. The curve is normalized so that 

M(r c ) = The error bars are: 

AAf(r) = ^(471)2 Jl /; Covin, riyfrtdndrz . (37) 

On the lower panel of Figure [7J we plot the potential of 
such an object as in Equation (l29l) alo ng with the top 
hat potential from lRudnick et al\ ( 20071 ) . 

Since the potential is only Equation (|29[) and in linear 
theory its time dependence is relatively simple $(r, r) = 

$(r)^^y, it is straightforward to integrate Equation \21\ 
numerically. In our calculations we placed the density 
fluctuation at the median redshift (z — 0.53) of the SDSS 
photometric LRG survey so th at we can compare o ur re- 
sult to actual measurements ()Granett et aLll2008D . In 
Figure |5] we plot AT versus the comoving radius (r c ) of 
the superstructure, we also plot AT according to Equa- 
tion (|2"5)l . Here, we used = 1. It is clear that Equation 
underestimates the ISW effect. 
In Figure O we used Equation (|2"6l to calculate <5 in . 
The PDF was measured in an SDSS LRG subsample at 
z = 0.28 median redshift (see 12.21 for details) and scaled 



Fig. 8. — AT for a photon traveling through the center of a su- 
pervoid at redshift 0.52 with Si n = 1 against its r adiu s is plotted 
(solid line). The dashed line comes from Equation J28J, an approx- 
imation using compensated profile. 



to the subsample described in IGranett et all {2008), 
which is located at z = 0.53. This means the scaling 
of the variance, S3 and S4 according to linear dynamics. 



The value for V. 



survey 



ca me fr om the properties of the 



survey in IGranett et all (|2008D . All three assumptions 
in Equations (|23H25|) about the density distribution give 
similar results. We also plot an estimate based on raw 
data without deconvolution. The dashed line is the ISW 
effect according to Equation (|2"51) . We plotted both the 
supervoids and superclusters. The intrinsic fluctuations 
of the matter density (see Figure [S]) and the uncertainty 
of the tail of the PDF add up. On the horizontal axis 
the scale is the linear size of the cell we measured CIC 



4. DISCUSSION 

We measure the hrst few moments of the matter PDF 
and give an estimate of the linear ISW effect owing to 
the largest density fluctuations on 100 /i _1 Mpc scales. 

First, we measure CIC in a volume- limited subsample 
of the SDSS DR7 spectroscopic LRG data. We com- 
pare the variance, S3 and S4 to their values in ACDM. 
Despite that we use the lowest order approximation we 
find agreement between data and prediction (Figure [S]) . 
However, it cannot be excluded that more complex mod- 
els can fit the data better especially on nonlinear scales. 
As can be seen from Figure [5l S3 and S4 differ slightly 
from the lowest order predictions possibly due to nonlin- 
earities, although this difference is not statistically sig- 
nificant. We approximate the tail of the matter PDF 
with analytic functions to get the density extrema. We 
calculate the expected radial profile of supervoids and su- 
perclusters with the condition that their average density 
inside a sphere is known. We estimate the average linear 
ISW signal of 50 of the most significant from each in a 
realistic survey. As can be seen from Figure |H1 it is plau- 
si ble that linear I SW c an produce the results presented 
in IGranett et a l. ( 2008). They used a compensating top 
hat filter with inner radius of 4° to get 7.9 ± 3.1{iK for 
clusters and — 11.3±3.1/Lt-ftf for voids. In comparison, our 





L(MPc/h) L(MPc/h) 



Fig. 9. — Expected ISW effect from the a vera ge of 50 supervoids 
LRG sample. The dashed line is Equation i28h On each panel we 

estimates for the temperature at the center of the same 
stacks from Figure are 5.5 ± 3.5 and —5.5 ± 3.5 in the 
case of a Gaussian PDF. We plot the linear size of a cube 
on the x-axis. The projection of a sphere with the same 
volume gives the corresponding angle. In case of L = 
150 /i _1 Mpc, this is 3°. 8. While our errors originate from 
the fluctuations of the ISW signal, theirs come from the 
primary CMB anisotropics. Thus they ar e independent 
and ou r calculation is fully consistent with lGranett et ali 
(2008). Our estimate is robust. It is not affected signifi- 
cantly by the details of the matter PDF. The error bars 
can be tightened if the volume of the survey is larger 
and if more images are stacked. The former would re- 
duce the cosmic error on the scales we study, while the 
latter would give a more accurate measurement of the 
average profile. 

We take one step toward cosmological parameter es- 
timation with calculating the expected linear ISW sig- 
nature of supervoids and superclusters. The next step 
can be to depart from the spherically symmetric model 
that we use for the sake of simplicity. Anisotropic fluc- 
tuations in the matter density around the center of a su- 
perstructure might give a quantitatively different error 
estimate. We also ignore any nonlinearities. We work 
with linear scales but we also probe the highest and low- 
est densities. The latter calls for a biasing model more 
complex than linear. We also ignore the nonlinear ISW, 



1 superclusters from a survey similar to the SDSS DR4 photometric 
a different approximation of the matter PDF. 

the Rees-Sciama effect (|Rees and Sciamalll968f ). It has 
been shown that i t is small compared to the linear part at 
low redshifts (see lCai et aZ.ll20d9l I2010T ). Another possi- 
ble improvement is to use general relativity (GR) instead 
of Newtonian. A model of compensate d voi ds based on 
GR is discussed in llnoue and Silk] (|2007| ) and llnoue et ali 
(2010). We also ignore the correlation between the ob- 
jects. 

Our result suggests that the void needed to produce 
detectable anomalies on the CMB is smalle r than pre- 
vious ly estimated. The CMB Cold Spot (jCruz et ali 
120051 ) has been considered consistent with a compen- 
sa ted void having S = —0 .3 and a radius of 200 /i _1 Mpc 
bv lSakai and Inouel (l2008h in ag r eemen t with the heuris- 
tic argmrient~^f^SdnIck TTalTlpOOl . From Figure M 
one can see that a void of similar size needs to have 
a much less significant underdensity in our Gaussian, 
non-compensated model. The signal from a top-hat 
potential is less than half of the signal of a realistic 
void. This makes the detecti on of such voids ha rder 
in today's sp a rse ca talogs (see iGranett et ali 120101 and 
iBremer et al\ (pOlOh ). 

We thank Ben Granett for his useful comments. The 
authors were supported by NASA grants NNX10AD53G 
and NNG06GE71G. 
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